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We consider a system consisting of two quantum dots connected by an s-wave superconductor in 
the presence of periodically varying electric or magnetic fields. The Floquet theory shows that there 
may be Floquet poor man's Majorana fermions (FPMMFs) in the high frequency region depending 
on the phase difference between the applied external fields to the two dots. Numerical results confirm 
this expectation and find in addition a lot of FPMMFs in the low frequency region. The FPMMFs 
survive for nonzero energy levels of the two dots and their interaction though the frequencies at 
which they emerge change and thus may be a promising candidate to be detected in experiments. 

PACS numbers: 74.45,+c, 85.35.Gv, 74.20.Mn 



I. INTRODUCTION 

Recent researches on Majorana's fermion whose anti- 
particle is itself bring exciting new physics in the fields 
of condensed matter physics and ultracold atoms es- 
pecially the Majorana fermion observed in the very re- 
cent experiments^ - — However, to definitely confirm the 
observed Majorana fermion, one still needs more careful 
analysis and other promising experiments, since factors 
such as disorder can also induce the experimentally ob- 
served zero-bias anormaly in differential conductance^ 
Therefore, many other promising designs for detecting 
the Majorana fermion have been proposed ] 10 ! 11 Among 
them, instead of quantum wires one can use quantum 
dot array to host Majorana bound state for the rea- 
son that the latter can overcome the disorder effects 
emciently i 12 i 13 Besides, easily-controlled quantum dots 
not only can be used for the detection of Majorana 
fermions in quantum wires, 14 but also may be helpful 
for non-Abelian statistics realizations. 15 

On the other hand, similar to Floquet topological 
insulators^ there is also a line to study the topo- 
logical superconductor system by periodically driving 
the chemical potential or external fields. The concept 
of Floquet Majorana fermions was introduced in the 
ultracold atoms ^ This line becomes more and more 
attractive^ - — since the rich physics may appear in the 
periodically changed systems such as nonzero energy Ma- 
jorana fermions in the Floquet space based on the fact 
that states can keep intact when Floquet quasienergies 
change by integral multiples of frequency (h is set to unit 
in this paper). Also, it has been pointed out that Floquet 
Majorana fermions follow the same statistics properties 
as ordinary ones. 21 

In this paper, we introduce Floquet poor man's Majo- 
rana fermions and study their properties on the basis of 
the very recently proposed interesting work of poor man's 
Majorana fermions (PMMFs). 24 The PMMFs share the 
same properties as the ordinary Majorana fermions, but 
lack topological protection. The PMMFs are so interest- 
ing because not only their detection is very feasible as 
there were already other experiments on the very similar 
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FIG. 1. (Color online) Schematic diagram of two quantum 
dots connected by an s-wave superconductor. The gates V g i 
and V g 2 tune the energy levels of the double dots. The mag- 
netic fields B± and Bi fix the spin directions of the electrons 
and provide nonzero values of A and A (see text). Two alter- 
nating circuits with local electrodes are connected to the dots 
to produce periodically varying chemical potentials. 



systems ^ 22 ! 23 but also they can locate at the double quan- 
tum dots separate without overlapping with each other. 

In the following, we first introduce the system we con- 
sidered in Sec. Inland recall the condition of the appear- 
ance of PMMF in Sec. HIl Then we describe in Sec. HVl 
how to use the Floquet theory to analyse the appearance 
of FPMMF when the adiabatic approximation is met. 
Sec. |Vl concentrates on numerical studies. We first cal- 
culate in Sec. IV Al the quasienergy and Floquet modes in 
the absence of interaction both in the low frequency and 
high frequency region in the Nambu space, even the adi- 
abatic approximation is violated. Then, we transfer from 
the Nambu space to the Fock Space to consider the inter- 
action effects in Sec. IV Bl The time evolution of Floquet 
modes is presented in Sec. V C followed by a summary 
in Sec. ED 



II. MODEL 

The system we considered is illustrated in Fig. HJ the 
double quantum dots are connected by an s-wave super- 
conductor. We assume for simplicity that there is only 
one state in each quantum dots, and the spin directions 
of the electrons are fixed by their respective magnetic 
fields Bi and B2. As a result, no double occupancy and 
intra-dot electron interactions exist and only the interac- 
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tion between the dots is considered. The Hamiltonian of 
the system is then, 



H = 6\d\di + e 2 d\d 2 + Ud\did\d 2 

+A (d\d 2 - + A (dJ4 - ^1^2 



(1) 



where d\ and are the creation and annihilation op- 
erators of electrons in dot i with the energy levels €i 
(i = 1,2), respectively, U is the interaction between the 
electrons, A is the co- tunneling hopping amplitude, and 
A represents the amplitude of cross Andreev reflections. 
The driving term produced by the electric or magnetic 
field is, 



H e (t)=f 1 {t) d\d 1 +f 2 (t) d\d 2 



(2) 



where fi = ^4cos(W) and i 2 — ^4cos(co>£ + uj is the 
driving frequency, t is the real time, A and $ are the 
driving amplitude and the phase difference. Note that 
the forms of the contributions of the electric field and 
magnetic field to the Hamiltonian are identical and the 
driven terms change in fact the chemical potential. How- 
ever, it is more easy to use the electric field in practical 
experiments as it is more easy to be controlled. 



III. POOR MAN'S MAJORANA FERMIONS 

We recall briefly in this section how the PMMFs ap- 
pear in such a system in Leijnse and Flensberg's paper. 24 
In their work, an inhomogeneous magnetic field are ap- 
plied to modify the values of A and A: A = Ao sin (0/2), 
A = Aocos(#/2), with being the angle between the 
the directions of magnetic fields applied on the two dots. 
Tuning the angle can then make the ratio A/ A equal unit. 
At this ratio, the solution of the Hamiltonian has two 
degenerated zero eigenvalues. The corresponding eigen- 
states are two well separated Major ana bound states: 
7i = (di + d\)/y/2 and 72 = i(d 2 — d) 2 )/\/2. 71 locates 
in dot 1 and 72 locates in dot 2. Thus tuning the ratio is 
one of the key points in forming Majorana bound states 
in the system. We shall show in the following how to 
tune the ratio using periodically varying external fields 
and reveal rich phenomena in this case. 



IV. ANALYTICAL RESULTS IN HIGH 
FREQUENCIES 

Let us introduce Floquet's theor y 25 ! 26 which is a the- 
ory to handle a system whose Hamiltonian is periodi- 
cally changed in time. According to this theory, for 
H(t + T) = H (t) with a period T = 2tt/uj, the wave 
function \ip(t)) of the Schrodinger equation, 



H\iP(t))=i^(t)), 




FIG. 2. (Color online) (a) Zero-order Bessel function and (b) 
its inverse as a function of frequency. In the high frequency 
region, a pair of FPMMFs appears (red circle) at the phase 
difference 2kir in (a), and the pair of FPMMFs disappears at 
the phase difference (2k + 1)ty in (b). We set A = 1 and the 
fixed ratio A/A = 0.8. 



can be expressed as 



\ip a (t)) = exp (-ie a t) \<p a (t)), 



(4) 



where e a is the quasienergy and \(p a (t)) is the periodic 
Floquet state satisfying \tp a (t + T)) = \(p a (t)). Substitut- 
ing this into Eq. (|3]), one obtains the Floquet equation, 



Hf\<p a (£)) = e a \(p a (£)) 



(5) 



with the Floquet Hamiltonian being defined as Hf = 
H — id/dt. One can also use the time evolution operator 



U(t +T,t ) = TexpH I t0 
Floquet equation as 



to+T 



H(t)dt) to express the 



U(t + T,t )\<p a (t )) = exp(-ie a T) \p a (t )>, (6) 

where T is the time order operator. Solving the Floquet 
equation instead of the original Schrodinger equation is 
the way of the Floquet theory to handle the nonequilib- 
rium problem. This is similar to the Bloch theory which 
handles a system with periodically changed Hamiltonian 
in space. Thus to diagonal the Floquet Hamiltonian is 
the main task for obtaining the quasienergy and the Flo- 
quet state. 

To this end, we use the rotating wave approximation, 27 

|{n^}, m) = exp ^imut — i J H e (t f )dt^ (7) 

to express the Floquet state, where to = 0, means 
a state in Fock space and the integer m is the quan- 
tum number of the so-called photon mode in the space. 
Absorbing or emitting a photon can be compensated by 
adding or removing a quasienergy equal to the photon en- 
ergy without changing the wave function of the original 
Schrodinger equation from Eq. (jl]). To ensure the orthog- 
onality of the Floquet states, one should define the inner 
product of the Floquet states over a period besides the 
traditional inner product, viz. 



1 f T 

(3) ((<Pa (t) \<Pj3 (*))) = t] o (<P« (t) \<Pj3 (t))dt = 5 at p. 



(8) 
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In the high frequency region, when A is relatively 
small and A/uj <C 1, the adiabatic approximation is 
met. 18 There is no mixture between different photon 
modes, since the system needs an energy uj to jump 
from one photon state to another photon state. Thus, 
one can only consider the contribution from zero pho- 
ton state, i.e., m = m' = 0, to calculate the matrix 
elements of the inner product of the Floquet Hamilto- 
nian, ({{ ni },m'\H f \{ ni },m)) = (({m}\H \{m})). This 
inner product is determined by the results of the com- 
mutation between H and H e . In our case, the possible 
nonzero commutation results from the hopping term and 
the cross Andreev reflection term, 



(d\d 2 - di4) , fi W 4 d i + f 2 (t) d\d 2 



= -XA 



(cos (out) — cos (ut + 1?)) (4^2 + ^4), (9) 



A (d\d\ - dida) , fi it) d\d x + f 2 {t) d\d 2 

= —AA (cos (ujt) + cos (uot + #)) (44 + • (10) 



Therefore, the phase difference d plays a key role in the 
commutation relation. If d = 2kir, for an integer k 
only the cross Andreev reflection term does not commu- 
tate with the driving term. The influence of the driving 
term then reflects on the matrix elements of the Flo- 
quet Hamiltonian through a renormalized cross Andreev 
coupling A — > AJo (2 A/a;), where Jo is the zero-order 
Bessel function*2!i2£ If d = (2k + l)7r, then the hop- 
ping term does not commutate with the driving term and 
one gets A — ^ \Jq(2A/uj). Thus, using the Nambu ba- 
sis = (di, d 2 , 4' 4)' the effective Floquet Hamiltonian 
can be expressed as 



with 
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(12) 



for $ = 2kn, and 
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(13) 



for <& = (2fc + 1)tt. 

The FPMMFs can emerge by tuning the frequency 
which changes the ratio A/ A. When ei = e 2 = 0, a 
properly tuned frequency makes the ratio A/ A = ±1; 
or, if we fix the value of A/A first, then, when the zero- 
order Bessel function equals the value of ±A/A, the so- 
lutions of hf have two degenerated zero energy solutions, 
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FIG. 3. (Color online) (a) Two quasienergies as a function 
of uj and (b) the four Floquet modes of the quasienergy (red 
curve in (a)) as a function of frequency at time zero: the 
orange curve, green curve, purple curve and blue curve rep- 
resent gi to #4, respectively, (c) The detailed behaviors of 
the Floquet modes in the low frequencies when the FPMMFs 
appear at time zero, (d) Details of quasieneries in the low fre- 
quency region, where the red circles represent the FPMMFs 
emerging. All with a phase difference 2/c7r. 
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FIG. 4. (Color online) Same with Fig. [3l but all at a phase 
difference (2k + 1)ty. Note that in (b) the orange curve is 
covered by the blue curve in the high frequency region. Both 
Fig. [3]and Fig. H are calculated in the Nambu space. 
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(1, 0, 1, 0) T and (0, -1, 0, 1) T for A/A or (0, 1, 0, 1) T and 
(— 1, 0, 1, 0) T for —A/ A. The operator expressions of the 
Floquet states are 71 = (d\ + d\)/y/2 and 72 = 2(^2 — 
4)/>/2 or 73 = (d 2 + 4)/V^ and 74 = z(di - 4)/V% 
respectively. However, in the high frequency region, one 
can only find one pair of FPMMFs, 71 and 72 for ex- 
ample, see the red circle in Fig. [2(a), if one fixed the 
value of the ratio to A/ A = 0.8 first. When the phase 
difference is (2k + l)7r, the inverse of zero-order Bessel 
function needs to cross the fixed value of the ratio for 
the FPMMFs to emerge. However, no solutions exist in 
this case. Therefore, for the same value of the ratio whose 
$ — 2/c7r, changing $ to (2fc + i)ir will extinguish the pair 
of FPMMFs from the quantum dots, see Fig. [2(b). 



V. NUMERICAL RESULTS 

In this section, we study the low frequency region that 
is more interesting, though one needs to include the ef- 
fects of the mixture of the different photon modes. We 
use the program QuTiP 28 to numerically calculate the 
quasienergy spectrum both in the high and low frequency 
region. The QuTiP finds the Floquet quasienergies and 
modes by calculating and diagonalizing the time evolu- 
tion operator U (to + T, to) m Eq. ©. 



A. No interaction 

We neglect the interaction U and set e\ = €2 = 0, 
A = 1, A = 0.3 and A = 0.24, then the ratio A/A = 0.8, 
and the numerical results can be compared with the anal- 
ysis in Sec.[lVl It turns out that in addition to confirming 
the above analysis, in the low frequency region many FP- 
MMFs appear under the condition of both $ = 2kir and 
$ — (2k + l)7r, even though it seems, at the first sight, 
that the FPMMFs would not emerge in the low frequen- 
cies for A/ A = 0.8 (see Fig. [2(a)). This indicates the 
contribution of the different photon modes. In Fig. [3(a) 
and Fig. [3(a), two different quasienergies cross over the 
zero quasienergy to constitute degenerate states. The 
Floquet modes at these points constitute the FPMMFs, 
see the zero time behaviors in Fig. [3(b), (c), and Fig. [4(c). 
We find that different ratios have different distributions 
of the FPMMFs in the low frequency region. However, 
in the high frequency region, when the ratio changes from 
0.8 to over 1, the only one pair of the FPMMFs disappear 
or appear corresponding to the phase difference chang- 
ing from 2kir to (2k + 1)tt or vice versa, see Fig. [3(b) 
and Fig. 2(b). If the energy level in any quantum dot is 
changed to a nonzero value, the FPMMFs still exist both 
in the high and low frequencies for $ = 2kir and in the 
low frequencies for d = (2k + 1)tt. Moreover, generally, 
there are no nonzero quasienergy FPMMFs in this case, 
since adding or reducing an energy uj will add or reduce 
a phase with exp (iut) to the Floquet modes and thus 
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FIG. 5. (Color online) (a) Two quasienergies as the function 
of uj and (b) the eight Floquet modes of the two quasienergies 
(red and black curves in (a)) as a function of frequency at 
time zero: the solid orange, green, purple and blue curves 
represent g± to g<± of the black color quasienergy, respectively. 
The dashed orange, green, purple and blue curves represent 
gi to #4 of the red color quasienergy, respectively, (c) The 
detail behaviors of Floquet modes in the low frequency when 
the FPMMFs emerging, (d) Details of quasieneries in the low 
frequency region, where the red circles represent the FPMMFs 
emerging. 
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FIG. 6. (Color online) Same with Fig. \E\ but all at a phase 
difference (2k + 1)tt. Note that in (b) the solid orange, blue 
curves and the dashed green, purple curves are overlapped in 
the high frequency region. 
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preventing the hermiticity of the Floquet modes from 
forming the FPMMFs. 

For the detail form of the FPMMFs, one can already 
see in Fig. [3^b) that #2(^0 (the green curve) and gi{uS) 
(the blue curve) of the FPMMFs in the high frequency 
region do not equal zero. So, the operator expression 
of the FPMMFs do not coincide with (1,0, 1,0) T in 
the above analysis. This means that there are some 
overlaps of the FPMMFs locating in different quantum 
dots, which are (g\,g2,9$,9±) with g^cj) g 3 (uj) and 
92^) = #4(^)7 respectively. Similar behavior shows in 
the FPMMFs in Fig. [3{c) and Fig.gfc). This is because 
the frequencies at which the FPMMFs emerge is not high 
enough. If we use a larger A and a smaller A to guaran- 
tee the high-frequency FPMMF in Fig. [3fb) appear in a 
higher frequency region to satisfy the adiabatic approx- 
imation, then the Floquet modes will be (#i, 0, #3, 0) T 
with gi = gs. In addition, we note that the only deviation 
of €1 away from zero does not change much the frequency 
of the appearance of the FPMMFs, but the deviation of 
62 from zero at the same time substantially lowers the 
frequency of appearance of the FPMMFs. Also, e\ 7^ 
or 62 7^ makes the overlap of the FPMMFs locate in 
different quantum dots. 



B. With interaction 

To consider the interaction effect between the double 
dots, we change the Hamiltonian from the Nambu space 
to the Fock space and use the basis {|00), |10), |01), |11)}, 
where \n1n2) means having m states in dot 1 and ri2 
states in dot 2, respectively. The positions of the four 
bases correspond to the Floquet modes {gi, #2, #3, #4} in 
the numerical process. The appearance of FPMMFs now 
dictates that the Floquet modes with different fermion 
parity (odd \a Q ) or even \/3 e }) have identical Floquet 
quasienergy. In other words, adding a fermion does not 
change the quasienergy of the system, which means the 
FPMMFs have zero quasienergies. The modes are ex- 
plicitly \a Q ) = (0, #2, #3, 0) T with nonzero #2 and gs and 
\Pe) — (gi, 0? 0? g±) T w ^h nonzero g\ and g±. We perform 
numerical calculations in the Fock space. We checked 
that the results in the Nambu space were reproduced 
well. We then let 61,62, and U all have nonzero val- 
ues and find that the FPMMFs still exist, see Fig. [5] at 
d = 2/ctt and Fig. [6] at 1? = (2k + 1)tt for ei = e 2 = 0.1, 
A = 1, A = 0.3, A = 0.24 and U = 0.1. In Fig. Efa), two 
Floquet quasienergies (red and black curves) have a cross 
point in the high frequency region and many cross points 
in the low frequency region. At the time zero, the corre- 
sponding Floquet modes are shown in Fig. [5fb): g\ (the 
solid orange curve) and #4 (the solid blue curve) of one 
quasienergy mode and g^ (the dashed green curve) and gs 
(the dashed purple curve) of the other quasienergy mode 
become zero while the others become nonzero starting 
from the high frequency cross point in Fig. [5ja). Thus 
there is a pair of FPMMFs appears in the cross point. 
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FIG. 7. (Color online) The Floquet modes as a function of 
the time within two periods. The frequencies of (a) and (b) 
are the values before and after, respectively, the appearance 
of the pair of FPMMFs in the high frequency region in Fig. [3] 
The frequencies of (c) and (d) are the values before and after, 
respectively, the appearance of the pair of FPMMFs in the 
high frequency region in Fig. [5] 



Similar situations happened when the other FPMMFs 
appear in the low frequencies as indicated by the red cir- 
cles in Fig. [5fd), see Fig. [5fc) and (d) for details. Note 
that not all cross points of the red and black curves cor- 
respond to the FPMMFs: one can see an example in (d) 
when the frequency is about 0.24. In a word, the interac- 
tion U also supports the existence of FPMMFs, though it 
lowers the frequency of the appearance of the FPMMFs 
in the high frequency and changes the distribution of the 
FPMMFs in the low frequency region. 



C. Time evolutions 

Note that the Floquet modes shown above are just 
at the time zero. To see the time evolution of the FP- 
MMFs, we perform numerical calculations using the Flo- 
quet modes at time zero as initial states. It turns out that 
the FPMMFs always exist in the whole time region, see 
Fig. [7| since the time evolutions of the Floquet modes 
before and after the frequency at which the FPMMFs 
appear guarantee that the Floquet modes at the nonzero 
time behaves similarly to the zero time. It is worth not- 
ing that when the adiabatic approximation is met in the 
higher frequency region, #2 (the green curve) and gs (the 
purple curve) in Fig. [7f a) will overlap and so will the g\ 
(the orange curve) and #4 (the blue curve) in Fig. Efb). 
However, different things happen with the solid orange 
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and the solid blue curves in Fig. [TJ^c) and the dashed 
orange and the dashed blue curves in Fig. [7(d). Since 
the nonzero energy levels and the interaction lowers the 
frequency at which the FPMMFs appear, the solid or- 
ange and the solid blue curves are hard to overlap and so 
are the dashed orange and the dashed blue curves. This 
means they enhance the overlap between the FPMMFs 
locating in different dots. 

VI. SUMMARY 

We have introduced the FPMMFs and studied their 
appearance in the Floquet space. A simple analysis has 
shown that there is a pair of FPMMFs emerging in the 
adiabatic approximation when the phase difference in the 
double dots is the same. Direct numerical calculations 



have confirmed this analysis and have found a lot of FP- 
MMFs at low frequencies owing to the contributions of 
different photon modes. Different to other Floquet Majo- 
rana fermions, there is no nonzero quasienergy FPMMFs. 
When the adiabatic approximation is not satisfied even 
in the high frequency region, there is overlap between 
the two FPMMFs appearing in the different dots. The 
nonzero energy levels of the two dots and their interaction 
do not prevent the emergence of FPMMFs. The differ- 
ent behaviors due to the driving phase difference between 
the double dots provide a method other than tuning the 
ratio between the hopping amplitude and cross Andreev 
reflection amplitude to observe Majorana fermions. The 
FPMMFs may thus be promising to be found in experi- 
ments. 

This work was supported by the CNSF (No. 10625420) 
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